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Abstract 

We discuss in detail the constraints on partial compositeness coming from flavour 
and CP violation in the leptonic sector. In a first part we present a formulation of par¬ 
tial compositeness in terms of a flavour symmetry group and a set of spurions, whose 
background values specify the symmetry breaking pattern. In such a framework we con¬ 
struct the complete set of dimension-six operators describing lepton flavour violation 
and CP violation. By exploiting the existing bounds, we derive limits on the composite¬ 
ness scale in different scenarios, characterised by increasing restrictions on the spurion 
properties. We confirm that in the most general case the compositeness scale should lie 
well-above 10 TeV. However, if in the composite sector mass parameters and Yukawa 
couplings are universal, such a bound can be significantly lowered, without necessarily 
reproducing the case of minimal flavour violation. The most sensitive processes are 
decays of charged leptons either of radiative type or into three charged leptons, ji ^ e 
conversion in nuclei and the electric dipole moment of the electron. In a second part 
we explicitly compute the Wilson coefficients of the relevant dimension-six operators in 
the so-called two-site model, embodying the symmetry breaking pattern discussed in 
our first part, and we compare the results with those of the general spurion analysis. 


1 Introduction 


If the solution to the gauge hierarchy problem is based on a new symmetry and not 
on antrophic considerations [1] or on special evolutions of the scalar sector in the early 
universe [2], new physics at the TeV scale is expected. In most of the existing models 
new degrees of freedom carrying flavour quantum numbers are present at the TeV scale, 
representing potential sources of flavour changing neutral currents (FCNC) and CP violation. 
So far both direct searches at the LHC and indirect searches in the context of precision tests 
and flavour physics have brought no conclusive evidence of new physics at the TeV scale. The 
negative outcome of the search for new physics in the flavour sector is particularly intriguing. 
Indeed a scale of new physics Ajvp as large as 10® TeV [3] is required by an effective operator 
analysis to preserve the good agreement between observations and theory predictions, unless 
the new flavour sector is highly non generic and involves specific mechanisms to suppress 
FCNC and CP violation at the desired level. 

The latter possibility is empirically supported by the huge hierarchies among fermion 
masses and fermion mixing angles, which can only be explained by some special dynamics. 
An effective mechanism suppressing FCNC and CP violation can be introduced by observing 
that in the electroweak theory the symmetry of the flavour sector is broken only by the 
Yukawa interactions. Minimal Flavour Violation (MFV) [4] is defined by the assumption 
that, even including new physics contributions, Yukawa couplings are the only source of such 
symmetry breaking. In MFV flavour effects from new physics are controlled and damped by 
the smallness of the fermion masses and mixing angles. In this framework data allow Ajvp 
to be considerably smaller, close to the TeV scale. MFV provides a useful benchmark for 
the discussion of the flavour sector, but it does not emerge as a unique framework from the 
known mechanisms aiming to explain the observed fermion spectrum [5] , or from the known 
models providing a solution to the gauge hierarchy problem. 

In this paper we reconsider the possibility that both the origin of fermion masses and the 
suppression of FCNC and CP violation are due to the mechanism of partial compositeness 
(PC) [6], as realized in the context of composite Higgs models (for reviews see [7,8]), and we 
perform a detailed analysis of flavour and CP violations in the leptonic sector. According to 
PC there are no direct couplings between the elementary fermions and the Higgs doublet. 
The Higgs doublet has potentially strong couplings to a composite sector, including, in the 
simplest case, a set of vector-like fermions with masses of the order of the compositeness 
scale. The SM fermions are mostly elementary and get their masses through mixing terms 
with operators of the composite sector, often modelled by vector-like fermions. 

An appealing realization of this idea involves anarchic Yukawa couplings in the composite 
sector. In this case the observed hierarchies between SM fermion masses and mixing angles 
are entirely due to the elementary-composite mixing terms. This is of particular interest, 
especially for the lepton sector, since the known pattern of neutrino masses and mixing 
angles as extracted from neutrino oscillation experiments [9] seems to support the idea of an 
underlying anarchic dynamics [10,11]. Higher dimensional operators describing low-energy 
FCNC and CP violations are depleted by both inverse powers of the compositeness scale 
and by the mixing terms, thus realising an efficient suppression mechanism known as RS- 
GIM [12]. Quantitative studies in concrete models show that in the anarchic scenario limits 
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from CP violation in the quark sector lead to Ajvp > 10 TeV [8], while the existing bound on 
the rate of /r —?• eq results in Ajvp > 25 TeV [8]. This strongly disfavors the anarchic scenario 
when the compositeness scale is of the order of 1 TeV. It is also known that PC at the TeV 
scale can satisfy the bound from flavour physics if Yukawa couplings in the composite sector 
are non-generic. For instance if we assume that such couplings are universal and, at the same 
time, that the only irreducible sources of flavour symmetry breaking are proportional to the 
SM Yukawa couplings, we reproduce exactly the MFV scheme [13] . 

In this paper we focus on the lepton sector. Several studies of lepton flavour violation 
(LFV) and CP violation in the lepton sector of composite Higgs models have been realised. 
Most of the analysis have been performed in the context of five-dimensional (5D) models 
with a warped space-time metric, weakly coupled duals to strongly coupled four-dimensional 
conformal theories, believed to provide a calculable framework for composite Higgs models. 
Explicit computations with anarchic Yukawa couplings have been carried out in refs. [14-20]. 
They get a lower bound on the masses of the first Kaluza-Klein modes of order 10 TeV. These 
bounds can be relaxed by requiring discrete [21-24] or continuous [25-27] flavour symmetries. 
The original motivation for discrete symmetries, tailored to approximately reproduce the tri- 
bimaximal mixing pattern in the lepton sector, has weaken after the precise measurement of 
the di 3 angle and the models considered in [21-24] require considerable corrections now. Con¬ 
tinuous non-abelian symmetries in the composite sector, broken by the elementary-composite 
mixing terms have been analysed in detail, especially in relation to FCNC and CP violation 
in the quark sector [28-32] . A review on flavour physics in 5D models with warped space-time 
metric can be found in ref. [33] . LFV and CP violation in the lepton sector have been inves¬ 
tigated also in SM extensions with extra vector-like heavy leptons [34,35], that can mimic 
the PC scenario, at least as far as the contributions from heavy fermions is concerned. 

In the present work we recast the framework of PC in terms of a generalized flavour 
symmetry and a suitable set of relevant spurions, much along the lines of refs. [12,28,32,36]. 
This involves a certain degree of model dependence, since both the flavour symmetry and the 
spurions are determined by the specific set of composite leptons, that we choose by following 
the criterium of minimality. In particular we work in the limit of vanishing neutrino masses, 
turning off the potential effects related to massive neutrinos. It is well-known that within 
PC, LFV and CP violation in the lepton sector are present also in the limit of massless 
neutrinos. By assuming that the set of adopted spurions are the only irreducible sources of 
flavour and CP violation, we can construct an exhaustive list of Wilson coefficients related 
to dimension six operators describing LFV and CP violating processes. At variance with the 
previous studies, we include for the first time all Wilson coefficients containing up to four 
powers of the spurions describing the elementary-composite mixing and we discuss their role 
in deriving the bounds on the compositeness scale. We also provide a complete list of the 
LFV Wilson coefficients that can be constructed in the limit of vanishing “wrong” Yukawa 
couplings. “Wrong” Yukawas in the composite sector are allowed by gauge symmetry, but 
do not contribute to SM lepton masses, at the leading order. They directly contribute to the 
dipole operators describing radiative decays of the charged leptons and setting such Yukawas 
to zero can relax the bounds on the compositeness scale. 

The general scope of our analysis is to check whether there are alternative solutions, 
beyond MFV, to reconcile PC at the TeV scale with the existing bounds on LFV and CP 


2 


violation. We also wish to verify if the anarchic scenario is completely ruled out or not. By 
exploiting the effective Lagrangian of our construction and the existing experimental bounds 
we estimate the limits on the new physics scale in several scenarios, where our spurious 
are subjected to a series of increasingly restrictive conditions. We confirm that in general 
the anarchic scenario is not compatible with new physics at the TeV scale and we provide 
examples of how PC can be realized at the TeV scale, without necessarily resorting to MFV. 

In the second part of our paper we consider as an explicit model realization of the flavour 
symmetry and its breaking pattern the so-called two-site model, first introduced in ref. [37]. 
Such realization contains explicitly vector-like leptons, implementing partial compositness in 
the lepton sector, as well as a set of spin-one resonances. By integrating out the states at the 
compositness scale, we evaluate the Wilson coefficients of the relevant LFV and CP-violating 
dimension six operators and we compare the results with those of the general spurion analysis. 

Our paper has the following plan. In section 2 we define the flavour symmetry and the 
set of spurious of our setup and we characterize the Wilson coefficients of the dimension-six 
operators relevant to LFV and CP violation in the lepton sector. In section 3 we perform 
a phenomenological analysis and we study the bounds on the new physics scale obtained by 
making different types of assumptions on the available spurious. In section 4 we recall the 
main aspects of the two-site model, that explicitly incorporates the features of the flavour 
symmetry breaking defined in section 2. In section 5 we collect our results on the dimension 
six operators obtained from the model by integrating out heavy fermions and heavy gauge 
vector bosons. In section 6 we present a phenomenological analysis of LFV in the two- 
site model. Finally we draw our conclusion. In the appendix B we show the result of our 
computation of the full one-loop contribution to the electromagnetic dipole operator in the 
two-site model. 


2 Effective field theory for lepton flavour violation 

As a first step, we have to choose the flavour symmetry group of our effective theory and 
its breaking terms. Throughout this paper we will work in the limit of massless neutrinos. 
The leptons are those of the SM, that is three copies of SU (2) doublets £ and singlets e. In 
MFV the flavour symmetry group of the leptonic sector is SU (3)^ x SU (3)e, corresponding 
to independent transformations made on £ and e. In our analysis we will instead assume 
a PC scenario. Charged leptons have no direct coupling to the Higgs doublet and acquire 
masses via mixing with vector-like heavy fermions. In this framework it is natural to assume 
as flavour symmetry group (focusing only on the non-Abelian part): 

Gf = SU{3f 

= SU{3)i X SU{3), X SU{3)l, x SU{3)l^ x SU{3)^^ x SU{3)^^ , (1) 

under which the lepton fields rotate in generation space in the following way: 

^Li iVe)ij£Lj , (Ve)ijeRj , ( 2 ) 

where Ve and W are elements of SU{3)i and SU{3)e, respectively. In other words, the SM 
leptons only transform under the SU (3)^ x SU (3)e component of the flavour group Gf, and are 
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invariant under the remaining SU (3)^ factor. Such a factor will be used to specify the spurious 
of our effective field theory. We introduce three sets of spurious. We need spurious (A, A) 
that mix SM leptons with heavy fermions. Moreover we allow for spurious (m, rh) describing 
the masses of the heavy fermions in the limit of unbroken electroweak symmetry. Finally 
the heavy fermion sector can interact with the Higgs doublet and this implies additional 
spurious {Y£,Y^). Notice that this set of spurious is the most general one compatible with 
our fiavour group Gf and with the assumption that the SM Higgs doublet directly couples 
only to the heavy sector. From these considerations the following transformation properties 
for our spurious can be deduced 


m ^ VL.mVl , A ^ ViAV^^ , 

m ^ , A ^ , 


^ Vl,Y*vI , 

Yl ^ , 


( 3 ) 


where it is evident from our notation which are the SU{3) factors of Gj that are involved in 
each transformation. An equivalent set of spurious is obtained by replacing the dimensionfnl 
quantities (A, A) with dimensionless combinations (X, A) = (Am“^, transforming 

as 

A ^ VeXVl , 

X VeXVl 

In such a case at the leading order (LO) SM charged lepton masses are generated by the 
operator 

= (iL'p)X Yi , ( 6 ) 

where (p is the Higgs electroweak doublet. 

It is well-know that, if all the spurious (A, X), (Yl,Y^) and (m, rh) are present at the same 
time and if the Yukawa couplings (Yl,Y^) are anarchic, a severe bound on the new physics 
scale applies. Indeed in concrete models belonging to the general class we are considering, 
the one-loop exchange of Higgs and heavy fermions leads to the following electromagnetic 
dipole operator: 

^JiLP)('y F)X Y^m-^Yl'm-^Yi X^e^ . (6) 

If such a contribution is present, and the Yukawa couplings (Yl,Y^) are assumed to be 
anarchic of 0(1), then the electromagnetic dipole operator and the mass operator are not 
aligned in fiavour space and the heavy fermion scale is bounded to be heavier than about 30 
TeV, to respect the bound on BR(/r —>■ ey). One way to eliminate this dangerous contribution, 
while maintaining non-vanishing lepton masses, is to assume Y^ = 0. We will come back 
to this assumption later in this section. For the moment we will adopt it as a working 
hypothesis. Our purpose is to analyse the fiavour violating contributions surviving in this 
limit and to estimate the corresponding bounds on the new physics scale. 

At variance with MFV, spurious with the dimension of a mass are present in our setup 
and some additional prescriptions are needed: 

• First of all we require that our operators are local in the spurious (X,X). These are 
mixing parameters that are generically treated as small and provide one set of expansion 
parameters for our spurion analysis. 


( 4 ) 
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• We also assume that the operators are local in the Yukawa couplings Y^, that we restrict 
in the range 1 < |Y^| < 47r. 

• Each power of occurs accompanied either by an Higgs electroweak doublet ip or by 
a factor l/dvr. 

• Masses of the composite sector are described by the spurious (m, fh) to which we add an 
additional parameter M, singlet under the flavour symmetry group, to describe masses 
of other composite particles, such as for instance a new set of vector boson resonances. 
These additional states are coupled to the Higgs doublet and to the heavy fermions 
with a strong coupling constant g*, in the range 1 < (/* < Att. 

• In our operators masses will always appear in negative powers, to allow for decoupling 
of all the operators with the scale of new physics. 

With this set of assumptions, the LO operator describing charged lepton masses is still Om 
in eq. (5). We stress that this operator provides the deflnition of the spurion Y^. Any 
polynomial of the type 




could replace Y^ in Om- With no loss of generality we can redefine as Y^ the particular 
combination of eq. (7) occurring in Om- 

The physical processes we are interested in concern lepton flavour violation in the charged 
lepton sector as well as the magnetic dipole moments and the electric dipole moments (EDM) 
of the charged leptons. To this purpose it is convenient to adopt an effective field theory 
description where the SM Lagrangian is extended by an appropriate set of gauge invariant 
operators depending on the SM fields [38]: 



( 8 ) 


where we have restricted our attention to the lowest dimensional operators relevant to the 
processes we are interested in, namely those of dimension six. Dots denote higher-dimensional 
operators. We list below a complete set of dimension six operators depending on lepton fields 
and on the scalar electroweak doublet ip [39]. We start with the dipole operators 


{Qew)ij = {iLiCr^''eRj)T^(pW^^ , 


( 9 ) 


and are the field strengths for the gauge vector bosons of SU{2) and U{1), respec¬ 
tively. The flavour structure of these two operators is the same and we will focus on the 
combinations 


{Qe'r)ij — COS 6w {QeB)ij “ Siu {QeW3)ij , 
{Qez)ij = Sin^vv {QeB)ij + COS dyi/ {QeW3)ij -, 


( 10 ) 
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where 6w is the weak mixing angle and QeW 2 , denotes the contribution to Q^w obtained 
setting to zero W^’^. The electromagnetic dipole operator is the only operator that gives 
a tree-level contribution to the radiative decays of charged leptons, when i ^ j. The diagonal 
elements, i = j, contribute to the anomalous magnetic moments and to the EDM of the 
charged leptons. We have operators bilinear in the Higgs doublet 

(<?)« = , 

, ( 11 ) 

, _ 

{Q<^e)ij iSP D^ Cjij^ 

After the breaking of the electroweak symmetry these operators modify the couplings of the 
Z boson to leptons, potentially violating both universality, through the diagonal terms, and 
lepton flavour, through the non-diagonal ones. There is a unique operator trilinear in the 
Higgs doublet 

which contributes, with a different weight, to both the masses and the Higgs couplings of the 
charged leptons. Finally we have four-lepton operators 

(Qu) ijmn 

{Qee)ijmn {^Rillx^Rj){^Rml'^^Rn) , (13) 

{Qle)ijmn Li'yLj^ Rm'y^^Rn) ) 

which can contribute to muon and tau decays into three charged leptons. We also consider 
dimension-six operators of the type llqq that can contribute to /r —?• e conversion in nuclei b 

~ {^Lilix^Lj){uLml'^ULn) 

^Qtq)i3 = dLn) 

i,Q£u,d)ij i,^Li'~ffi^Lj^i,^Rm'y^^Rn) i iJ,^Lj)i,djijji'y^dji^ , 

iQeq)ij {^Ri'yiJ,^Rj'){QLm'y^QLn) Q {'^id'), (13) 

i,Qeu,d)ij i,('Ri'yfi('Rj^i,^Rm'y^^Rn) i fj.^Rj^idRm'y^dji^ . 

There are other 3 independent operators of this type [39] , but the chosen subset is sufficiently 
general for the purposes of the present discussion. Notice that each operator carries flavour 
indices. Hermiticity of the effective Lagrangian is guaranteed either by appropriate symmetry 
properties of the Wilson coefficients under transposition of the family indices or by addition 
of the hermitian conjugate operator. 

Our aim is to estimate the Wilson coefficients of these operators, by expressing them 
in terms of the spurions using the set of rules described above. We expand each Wilson 
coefficient in powers of the mixings {X,X) and the Yukawa coupling Y^. Since the spurions 
{X,X) control lepton masses, they are expected to be small, of order (O.l/Y^) at most. We 

^The operators and are linear combinations of the operators and of 

ref. [39]. 
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Ny 

Nx 

HF 

HB 

1 

2 

XY*X^ 

X{cc<)-'YiX^ 

1 

4 

XY*X^XX^ 

XX^XY^X^ 

X{cc^)-^Y*X^XX^ 

XY^{cN)-^X^XX^ 

XY^X^X{c^c)-^X^ 

X{cc^)-^X^XY*X^ 

XX^X{cc^)-^Y*X^ 

XX^XY^{cN)-^X^ 

3 

2 

XY*Y*^Y*X^ 

X{cc^)-^Y*Y*^Y*X^ 

YY^(ctc)-W^^Y^Yt 

XYp^*\cc^)-^Y*X^ 

XY*Y*^Y*{c^c)-^X^ 

3 

4 

XY^X^XY^^Y^X^ 
XY*Y*^ X^ XY^X^ 
XYp^*^Y*X^XX^ 
XX^XY^Y^^Y^X^ 

X{cc^)-^Y^X^XYpY^X^ 

Y(cct)-W^Y^^YtYY^Yt 


Table 1: Spurion combination Ciji, in a matrix notation, for the lepton bilinear luiCi^ijCRj. 
Ny and Nx are the orders of the expansion in and {X,X), respectively. We restrict the 
list to Ny < 3 and Ny < 4. For convenience we distinguish spurion combinations depending 
on composite fermion matrices c and c (column HB), from combinations not involving c and 
c (column HF). 


will stop the expansion in {X,X) at the fourth order. In the expansion we will go up to the 
third order in Y^, since in the anarchic scenario trilinear combinations of Yukawa couplings are 
in general misaligned with respect to linear ones. Higher orders in Yukawa couplings do not 
bring any new qualitative feature in our analysis. Finally the correct dimension is provided 
by negative powers of M or {m,fh). Since the operators under study have dimension six, in 
practice we have two classes of operators, those suppressed by 1/M^ and those suppressed 
by the heavy fermion masses. Each of these two classes refers to a specific decoupling limit. 
When M ^ |m|, |m|, the heavy bosons decouple first and the operators are suppressed by 
negative powers of {m,m). We call this heavy boson (HB) case. In the opposite limit, 
M -C |m|, |m|, we have a fast heavy fermion decoupling and the operators are suppressed 
by the smaller scale M. We call this heavy fermion (HF) case. When the two scales M 
and (m, m) are comparable, the Wilson coefficient can be a generic function of the ratio of 
the two scales. For the present discussion the two limiting cases are sufficient to capture 
the behaviour of the system. For our spurion analysis it is convenient to rewrite the mass 
matrices m, in in this way 


m = mo c , rh = rriQ c , (15) 

where mo is a flavour-independent mass parameter, while the flavour dependence is carried 
by the dimensionless matrices c, c. To facilitate the identification of the relevant Wilson 
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Ny 

Nx 

HF 

HB 

0 

0 

1 

1 

0 

2 

xxt 

X{cc^)-^X^ 

0 

4 


X{cc^)-^X^XX^ 

XX^X{cc^)-^X^ 

2 

2 

XY*Y*^X^ 

X{cc^)-^Y*Y*^X^ 

XY*{c^c)-^Y*^X^ 

XY*Y*\cc^)-^X^ 

2 

4 

XY*X^XYpX^ 

XY^Y^^X^XX^ 

XX^XY^Y^^X^ 

X{cc^)-^Y*X^XYpX^ 

X{cc^)-^Y^Y^^X^XX^ 


Table 2: Spurion combination Cij^^ in a matrix notation, for the lepton bilinear ^Li{CLL)ij^ Lj- 
Ny and Nx are the orders of the expansion in and {X,X), respectively. We restrict the 
list to Ny < 2 and Ny < 4. For convenience we distinguish spurion combinations depending 
on composite fermion matrices c and c (column HB), from combinations not involving c and 
c (column HF). 


coelRcients, it is useful to identify the combinations of spurious that fit lepton bilinears 

^Li{Clji)ijejij , ^Li{ClL)ijlLj 1 ^Ri{Cjiji)ijeRj . (16) 

In tables 1, 2 and 3 we collect such combinations, constructed with the rules outlined above. 


Ny 

Nx 

HF 

HB 

0 

0 

1 

1 

0 

2 

xxi 

X{c^c)-^X^ 

0 

4 

XXtXXl 

X{c^c)-^X^XX^ 

XX^X(c^c)-^X^ 

2 

2 

XYpY*X^ 

X{c'^c)-^Y*^Y*X^ 

XY*^[cc^)-^Y*X^ 

XY^^y^(ctc)-ixt 

2 

4 

XY*^X^XY*X^ 

XY*^Y*X^XX^ 

XX^XY^^Y^X^ 

X{c^c)-^Y*^X^XY*X^ 

X{c^c)-^Y^^Y^X^XX^ 


Table 3: Spurion combination in a matrix notation, for the lepton bilinear eRi{Cfi^ijeRj. 
Ny and Ny are the orders of the expansion in and {X,X), respectively. We restrict the 
list to Ny < 2 and Nx < 4. For convenience we distinguish spurion combinations depending 
on composite fermion matrices c and c (column HB), from combinations not involving c and 
c (column HF). 
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2.1 Dipole operators 

The dipole operators in eq. (9) involve the lepton bilinear luiCi^ijCRj and their pos¬ 
sible Wilson coefficients are the {CiR)ij combinations fisted in table 1. We expect that the 
dipole operators are loop-generated in pertnrbation theory. The naive loop snppression factor 
l/(167r^) is not present in table 1 and should be included in the Wilson coefficient 

iCe-f,z)ij = ■ ( 17 ) 

The new physics scale A associated to dipole operators is M in the HF case (M <C \m\, \m\) 
and mo in the HB one (M ^ \m\, \rh\). 

The first important outcome of our analysis is that there is a large set of potentially 
lepton flavour violating combinations, beyond that of eq. (6). The actual appearance of 
these combinations in concrete models containing our set of spurious will depend on the 
specific dynamics of the model under consideration. We will discuss the phenomenological 
implications of the new structures for the Wilson coefficients {Ce^^z)ij in section 3. 

2.2 Scalar operator 

The scalar operator {Qe,f)ij = has exactly the same flavour structure of the 

dipole operators and we can directly read from table 1 the list of possible Wilson coefficients, 
provided we have at least three Yukawa couplings, Ny > 3, since this operator contains three 
Higgs doublets. Up to an overall flavour-independent coefficient, which is expected to be of 
order one, we have: 

= iClR)iJ (iVy>3) . (18) 

Notice that there is no loop suppression in this case, since this operator can be generated at 
tree level. As for the dipole operators, the scale of new physics A is M in the HF case and 
mo in the HB one. 

This operator, together with the mass operator of eq. (5), contributes to both masses and 
Yukawa couplings of charged leptons. If the Wilson coefficients are not exactly aligned in 
flavour space, we have flavour violating decays of the Higgs. Even in the case where there is 
a perfect alignment among all Wilson coefficients, the overall strength of Yukawa couplings 
is altered compared to the case of the SM, where it is completely fixed by the fermion mass 
and by the electroweak VEV. We will discuss the impact of these modifications in section 3. 

2.3 Vector operators 

In the operators {Q^e)ij of eq. (11) new flavour structures arise. The corre¬ 
sponding Wilson coefficients can be read from table 2, for and table 3, for (Q,pe)ij, 

but excluding the case Ny = 0, since the vector operators are bilinear in the Higgs doublets 
and this requires at least two powers of the Yukawa couplings. We have 

(^(^f )b ^ ; iCipe)ij = {Crr)^ , {Ny > 2) . (19) 

There is no loop suppression, in general, and the new physics scale is identified as for the 
previous operators. We see that, in general, the Wilson coefficients for the vector operators 
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are not diagonal in the mass basis. This leads to violation of flavour in the couplings of the 
Z boson, with consequences that we discuss in section 3. 


2.4 Contact operators 

Finally we consider the contact operators (QiOpmn, iQee)ijmn and {Qie)ijmn given in 
eq. (13). In these operators we recognise combinations of the flavour structures already 
discussed. We expect the following possible factorisations: 


(CO 


ijmn 


(C-ll) in ^ {CLL)mj 


( 20 ) 


(Ce) ijmn 


i^RR)ij {CRR)mn 

{^RR)in ^ {^RR)mj 

{CLL)ij {CRR)mn 
{ClR)in X {Cljjmj 


( 21 ) 

( 22 ) 


Since the contact terms contain no Higgs doublets, we have no restrictions on Ny. For 
the same reason, each bilinear in the Yukawa coupling should be accompanied by a loop 
suppression factor l/(167r^) Concerning the new physics scale A, when the two factors C 
come from the HF column we have A = M. When one the factors come from the HF column 
and the other from the HB column, we have A = mo. Similarly, for the llqq operators we 
have: 




iq J^J — y^LLJij 

{Ceq)ij = {CRR)ij 


{Ciu4)ij — {CRR)ij , 
{Qeu,d)ij (Crr) ij 


(23) 

(24) 


with no restrictions on Ny. 


2.5 Stability of the solution = 0 

Since we are interested in the scenario where is nearly vanishing, a legitimate question 
is whether and to which extent such a limit is stable under quantum corrections. To this 
purpose it is better to distinguish the two regimes v <C mo M and u <C M <C mo. When 
V <C mo <C M, we can consider the following spur ion combinations 


ml 

M2 



ml 

M2 


^ ^R 


(25) 


that behave as effective Yukawas of type ¥[ since they have the same transformation prop¬ 
erties as ¥[. We expect that terms like those in (25) arise, in perturbation theory, through 
threshold corrections induced by the one loop exchange of heavy gauge bosons and heavy 
fermions, from which we estimate 


Y* 


9* < Jy*c^ 

167r2 M2 ^ 


V* 


h 9* mp Ay*y*^V* 

167r2 M2 « * * 


(26) 


2lf contact terms originate from dim-8 operators, the loop factor can be effectively replaced by jbZ. 
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where is a coupling constant of the strong sector and k is either an additional loop factor 
stemming from a Higgs loop, or a factor coming from the electroweak VEV, if 

the two Yukawas are attached to external Higgs legs. The contributions of eq. (26) are small 
if there is an hierarchy between heavy fermions and heavy gauge boson masses. They are 
also suppressed in the semi-perturbative regime g*,Y^ 1. Similarly, exchange of ordinary 

gauge bosons and heavy fermions lead to effective Y£ couplings of the type 


Y* 


9^ '^w,z 

IGtt^ ttIq 


-1 


Y, 


s-i 


R'- 


Y* 


9^ ^w,z 

IGtt^ ttIq 



(27) 


These contributions are naturally suppressed by the ratio can be easily kept 

at the percent level. 

In the other regime, u <C M <C mo, beyond the combinations of eq. (27), we can also 
consider 


m2 


rrin 


'YrC 


m2 


mn 


-1 


ySYr 




YrC 


-1 


( 28 ) 


which also transform as the spurion Y^. These combinations formally decouple in the limit 
of infinitely large mo. There are however also combinations that do not decouple, such as 


c 


-1 




(29) 


We have checked, through a one-loop computation in the two-site model that we will con¬ 
sider in section 4, that this kind of contributions are generated. They cannot be paramet¬ 
rically suppressed by mass ratios and the corresponding Y^ can be depleted only in the 
semi-perturbative regime p*, Y^ ~ 1. 

Notice that for mo ~ M, the effective Yukawa couplings arising from eqs. (26,28), are 
insensitive to the overall mass scale. In summary, in the regime of heavy gauge bosons much 
heavier than heavy fermions, the effective Y^ can remain close to the percent level even in a 
strongly coupled regime, while in the opposite regime the solution Y^ ~ 0 is typically stable 
under quantum corrections (up to the % level) only if we assume the semi-perturbative regime 
9 *-: Yr ~ 1- 


3 Model independent bounds on the scale of new physics 

Model-independent studies of LFV processes have already been done in the literature 
[53-55]. In this section we collect the bounds on the Wilson coefficients of the Lagrangian of 
eq. (8) and discuss their impact on our spurion analysis. The main bounds on the coefficients 
of the dipole operators and O^z are given in table 5, showing the results obtained in 
ref. [54, 55] and derived from ref. [56] in the case of /i“Au ^ e“Au. They come from 
the present limits reported in table 4. The amplitudes for these processes get a tree-level 
contribution from the off-diagonal elements of the electromagnetic dipole operator while 
contributes at one loop. Each bound has been derived by assuming a single non-vanishing 
Wilson coefficient at the time. This also applies to all the bounds discussed in this section. 
Bounds on the coefficients C'g*^ are equal to the bounds on the coefficients The 
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LFV Process 

Present Bound 

Future Sensitivity 

/i eq 

5.7 X 10-1^ [40] 

^6x 10-1^ [41] 

fjj —y 3c 

1.0 X 10-12 [42] 

10-1® [43] 

/i“Au ^ e“Au 

7.0 X 10-1^ [44] 

? 

fx~Ti — )■ e“Ti 

4.3 X 10-12 [45] 

? 

H A1 — y e A1 

— 

10-1® [46,47] 

r ^ eq 

3.3 X 10-® [48] 

~ 10-® - 10-9 [49] 

T ^ 

4.4 X 10-® [48] 

- 10-® - 10-9 [49] 

T —7- 3e 

2.7 X 10-® [50] 

- 10-9 - 10-1° [49] 

T —3p 

2.1 X 10-® [50] 

~ 10-9 - 10-1° [49] 

Lepton EDM 

Present Bound 

Future Sensitivity 

de{e cm) 

8.7 X 10-29 [51] 

? 

d^(e cm) 

1.9 X 10-19 [52] 

? 


Table 4: Present and future experimental sensitivities for relevant low-energy observables. 



|C| (A 


: 1 TeV) 

A (TeV) (|C| = 1) 

LFV Process 


2.5 

X 

10-10 

6.3 X 10^ 

/i —7- eq 

^67 

4.0 

X 

10-9 

1.6 X 10^ 

jjj —y 3c 

^67 

5.2 

X 

10-9 

1.4 X 10^ 

fj,~Au —7- e-Au 

/^re 

^67 

2.4 

X 

10-® 

6.5 X 102 

r ^ ej 


2.7 

X 

10-® 

6.1 X 102 

T —>■ /rq 


1.4 

X 

10-" 

2.7 X 103 

fi ^ ej [l-loop] 

/^re 

^eZ 

1.3 

X 

10-^ 

28 

T —7- eq [l-loop] 

C7z 

1.5 

X 

10-3 

26 

r ^ pq [l-loop] 


Table 5: Bounds on off-diagonal Wilson coefficients and ^ from ref. [54,55]. The 

bounds from /r“Au —>■ e“Au have been derived from ref. [56]. In the second column we list 
the upper bound on ^1 assuming A = 1 TeV, while in the third column we fix ^1 = 1 
and we list the corresponding lower bound on A, in TeV. Bounds on the coefficients 6*^*2 
are equal to the bounds on the coefficients 2 - 


diagonal elements of the dipole operators contribute to electric and magnetic dipole moments 
of the charged leptons. From the present bounds reported in table 4 we have [54] 

< 3.9 X 10-^2 ^ < 8.4 x lO”" . (30) 

Given the current deviation Aa^ = in the muon anomalous magnetic moment 

a, = {g-2)j2 [57,58] 

Aa^ = (29 ± 9) X 10“^° , (31) 

we would need ^ 

fle(C«‘)(ii^) =1,2x10-= (32) 
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\C\ (A = 1 TeV) 

A (TeV) ( |C| = 1) 

LFV Process 

cue 

^ e^p 

8.4 X 10-^ 

109 

fi ^ ej [2-loop] 


0.33 

1.7 

T —>■ ey [2-loop] 

( 7 ta » 

^ CLP 

0.37 

1.6 

T ^ jj'j [2-loop] 


Table 6: Bounds on off-diagonal Wilson coefficients {Ce^)ij from ref. [59-61]. In the second 
column we list the upper bound on assuming A = 1 TeV, while in the third column 
we fix I(7*^ I = 1 and we list the corresponding lower bound on A, in TeV. Bounds on the 
coefficients (7^* are equal to the bounds on the coefficients (7*^. 


to account for the central value of the discrepancy. 

Also the scalar operator O^ip is mostly bounded by the limits on radiative lepton decays 
[54,55,59-61]. The scalar operator contributes to lepton masses and to higgs couplings with 
a different weight: 

Cy = M-ij &Li^Rj + Tp h eu^Rj + h.c. , (33) 

where 


Mij = 


SM 


-Vij + 


2A2 


{Cep) 


'f'J 


J ^/2 ’ 


v/2 


SM 


-Vij + 


3^2 


{Cep) 




(34) 


and are the Standard Model Yukawa couplings. Radiative charged lepton decays con¬ 
strain the off-diagonal elements of Tp in the basis where the mass matrix Adp is diagonal. 
To convert these constraints in bounds on the Wilson coefficients {Cep)ij it is convenient to 
work in the basis where the SM couplings are diagonal and expand the unitary matrices 
that diagonalise M. in powers of j IS?. We found that, to first order in this parameter, the 
off-diagonal elements Tp {i j) in the lepton mass basis are given by: 


= 


y2A' 


iC 


ecpjij 


(> ^ j) 


(35) 


By translating the bounds on Tp {i ? j) given in ref. [59-61] into bounds on {Cep)ij^ we 
get the results shown in table 6. The bounds on the coefficients {Cep)ji are equal those on 
the coefficients {Cep)ij- These bounds are dominated by two-loop contributions of the corre¬ 
sponding operator to the radiative lepton decay, through Barr-Zee type diagrams, assuming 
that the top Yukawa coupling is as in the SM. One-loop contributions to charged lepton 
radiative decays and tree-level contributions to £ 3£' decays lead to less severe bounds 

than the ones given in table 6. 

Coming to the vector operators {Qpf^)ij, {Qpe)ij, they lead to lepton flavour violating Z 
decays, but the corresponding limits on the Wilson coefficients, assuming A = 1 TeV, are of 
order 10% [54]. Through one-loop diagrams they also contribute to radiative decays of the 
charged leptons [54,55]. It turns out that the most restrictive bounds come from the processes 
/i“Au —7- e“Au and £ —?• 3£' whose branching ratios satisfy the experimental limits of table 4. 
We collect the corresponding bounds in table 7. Also the contact operators can contribute 
to /r“Au —>■ e“Au, £ —>■ 3£' and, through one-loop diagrams, to the radiative decays of the 
charged leptons. The most significant bounds are given in table 8. 


13 

















|C| (A 


1 TeV) 

A (TeV) (|C| = 1) 

LFV Process 

VTVTsyr 

Xipl )ne 

3.7 

X 

10-® 

164 

fjj —y 3c 

Xipl Ve 

5.0 

X 

10-6 

447 

li~Au —>■ e-Au 

Jre 

1.5 

X 

10-2 

8.3 

r ^ 3e 

vSw 

1.2 

X 

10-2 

9.0 

T — y Sjjj 

Cfie 

(fie 

3.9 

X 

10-6 

160 

jjj —y 3c 

cue 

^ Lpe 

5.0 

X 

10-6 

447 

fj,~Au — 7- e-Au 

r^re 

1.5 

X 

10-2 

8.1 

r ^ 3e 

^ Lpe 

1.3 

X 

10-2 

00 

bo 

T —y 3/i 


Table 7: Bounds on off-diagonal Wilson coefficients and from ref. [54,55]. The 

bounds from /r“Au —>■ e“Au have been derived from ref. [56]. In the second column we list 
the upper bound on the Wilson coefficients assuming A = 1 TeV, while in the third column 
we set to unity the coefficients and we list the corresponding lower bound on A, in TeV. 



\C\ (A = 1 TeV) 

A (TeV) (Id = 1) 

LFV Process 

^llee 

2.3 X 10-6 

207 

fjj —y 3c 

/^eree 

^llee 

9.2 X 10-6 

10.4 

T —>■ 3e 

^ll,ee 

7.8 X 10-6 

11.3 

r ^ 3/r 

^/j.eee,ee/j,e 

3.3 X 10-6 

174 

/jj —y 3c 

^/j,/j,/j,e,e/j,/j,/j, 

2.1 X 10-^ 

69 

/r —>■ ey [l-loop] 

^IJ,rre,err/j, 

1.2 X 10-6 

289 

/r —>■ ey [l-loop] 

^eree,eeeT 

^le 

1.3 X 10-2 

00 

bo 

T —>■ 3e 


1.1 X 10-2 

9.5 

T — y 3/i 

^ in 

2.0 X 10-6 

707 

fjj All —y c All 


1.8 X 10-6 

745 

fi~Au —)■ e-Au 


9.2 X 10-^ 

1.0 X 106 

fjj All —y c All 

^iu,eu 

2.0 X 10-6 

707 

fi~Au —)■ e-Au 


1.8 X 10-6 

745 

fjj All —y c All 


Table 8: Bounds on coefficients from ref. [54]. The bounds from /i“Au e“Au have 

been derived from ref. [56]. In the second column we list the upper bound on the Wilson 
coefficients assuming A = 1 TeV, while in the third column we set to unity the coefficients 
and we list the corresponding lower bound on A, in TeV. 


We can translate the bounds collected in tables 5-8 into limits on the masses M and ttiq 
of our spurion analysis. To do this we should determine or make some assumptions on the 
parameters c, c, X, X, and ¥£, through which we can express all the Wilson coefficients, 
as explained in section 2. By exploiting the flavour symmetry of our setup, we see that it is 
not restrictive to work in the basis where the mixing matrices X and X are diagonal, real and 
non-negative and we will adopt this choice, unless otherwise stated. The LO mass matrix of 
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the charged fermions 


V 


(36) 


m, = XY^ 



5 


is diagonalized by a bi-unitary transformation 


(L mi R^).. = mAj . (37) 

Assuming that is anarchic, it is straightforward to find that 


/ 1 

X 1 /X 2 

X 1 /X 3 \ 

( .K 

X 1 /X 2 

X 1 /X 3 \ 


X 1 /X 2 

1 

X 2 /X 3 , 

X 1 /X 2 

1 

X 2 /X 3 

, (38) 

\ X 1 /A 3 

X 2 /X 3 

1 J 

\ X 1 /A 3 

X 2 /X 3 

1 ) 



where factors of order one have been omitted from the matrix elements. Similarly the lepton 
masses are approximately given by 


mi ^ Xi{Y^)iiXi^ . (39) 

The most favorable scenario to minimize FCNC effects is realized when Xi = Xi, which we 
assume, for the time being. We have 


Xi = w = 


\/2i 


rui 




1/2 


(40) 


which will be used in our estimates. 

If c, c, and are all generic matrices, the dipole operator, eq. ( 6 ), leads to the 
well-known limit 

mo ^ > 33 TeV, (41) 

where (Y) denotes a suitable average of the Yl j^ couplings and similarly for (c). Such a 
bound can be evaded only if special features of the fiavour parameters c, c, X, X, Y^ and 
Yl are adopted. One possibility is to recover the framework of minimal flavour violation 
(MFV), as proposed in ref. [13]. This can be done by assuming c, c, X, X, Y^ and Y^ 
all proportional to the identity matrix, except either X (right-handed compositeness) or 
X (left-handed compositeness). In this case, when neutrino masses are neglected, the only 
spurion that breaks the flavour symmetry becomes proportional to the charged lepton Yukawa 
couplings, exactly as in MFV. A choice of basis where such a spurion is diagonal is always 
possible and LFV is only present when neutrino masses are turned on. In the latter case one 
can reconcile LFV with a scale of new physics close to the TeV scale, however we think it is 
interesting to explore other options allowing for a TeV scale A. 

We start by taking Yl = 0. To guarantee the stability of this condition we should 
also assume either a perturbative regime, where |Y^| < 1, or the hierarchy M !$> mo, as 
explained in section 2.5. In this case the dipole operator is dominated by the contribution 
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(Ny,Nx) = ( 1 , 2 ) (see table 1 ) 


r xr^(5tg)-ixt 

^67 ^ I ml 

Is? I X(ccl)-^Y^Xl 

I ml 

and the bound on /i —?• 67 leads to mQ{c) > 33 TeV. As a next step we consider the case of 

= 0 and universal heavy fermion masses, namely c and c proportional to the unit matrix. 
With this assumption the Wilson coefficients (Ay, Ax) = ( 1 , 2 ) of the dipole operator are 
diagonal in the mass basis, to LO. This however is not the case, in general, for the coefficients 
(Ay, Ax) = (3, 2), which are comparable in size to those of eq. ( 6 ). Therefore, even working 
in the limit of vanishing “wrong” Yukawa coupling ¥£, and universal heavy fermion masses, 
we expect that with anarchic Y^ the bound of eq. (41) still applies. To allow for a lower new 
physics scale A, we are lead to postulate that also the Yukawa couplings Y^ are universal. 
More precisely we define a new scenario, that we call intermediate flavour violation (IFV), 
by the following assumption: 

IFV scenario: 

The Yukawa couplings ¥£ are negligible and a ehoiee of basis exists where c, c and Y^ 
are simultaneously proportional to the identity matrix. 

From the discussion in section 2.5, we know that a non-vanishing Yf can be generated. 
Since such a Yf will only depend on c, c and Y^, it will be universal too. Admitting a 
universal Yf would not change our conclusions and, for simplicity, we keep the assumption 
Yf = 0. Moreover, from now on we set to unity the matrices c and c, absorbing their effect 
in the overall scale tuq of the heavy fermion masses. The remaining coupling Y^ is described 
by a single parameter y: 

Yi = yt . (43) 

At variance with MFV, in IFV the sources of flavour breaking are both X and X. In the 
basis defining IFV the mixing matrices X and X are not diagonal, in general. They are 
generic complex matrices. At the LO, the lepton mass matrix is proportional to the product 
XXT When we move to the charged lepton mass basis, in general, only the product XX^ 
becomes diagonal, not X and X individually and they both can lead to LFV. 

By exploiting the freedom related to our flavour symmetry, we can provide an alternative, 
but equivalent, description of the IFV scenario. We can choose a basis where X and X are 
diagonal, real and non-negative. In this case by means of the symmetry transformations of 
eqs. (2-4) we can still maintain c and c proportional to the identity, but the matrix Y^ becomes 
a generic unitary matrix. In such a basis the lepton mass matrix XY^X^ is non-diagonal 
and LFV is now ascribed to the interplay between Y^ and X, X. 

Within IFV the special case X = X, which was previously assumed to minimize FCNC 
effects, forbids any LFV effect. Indeed, if A = X, there exist a basis where both X and the 
mass matrix AY^Af are diagonal at the same time. As long as neutrino masses are neglected 
there is no source of LFV and all the Wilson coefficients are diagonal in flavour space. The 
only difference with respect to MFV is that the spurion A is not proportional to the charged 
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lepton Yukawa couplings, but to their square root. Therefore in discussing IFV we consider 
the general case where X and X are not equal. 

In IFV no LFV is generated from the Wilson coefficients of the dipole and scalar operators 
with Nx = 2, to LO. These coefficients are of the type XAX^, where A is any combination 
of c, c and and they are automatically aligned to the charged lepton mass matrix, to LO. 
To estimate the other Wilson coefficients we should specify the choice of X and X. Rather 
than scanning for the most general possibility, here we provide an example. We assume that, 
in the basis where X and X are diagonal, their elements are given, up to coefficients of order 
one, by 


Y = —diag(A^A^A3) 

y/y 


Y = —diag(A^A^l) 

y/y 


(44) 


where A ~ 0.22 is a small parameter of the order of the Cabibbo angle. The choice of elements 
of nearly the same order of magnitude for X is motivated by a possible role that the matrix 
X can play in describing large lepton mixing angles, once neutrino masses are turned on [5]. 
The value A^ ^ A^ is chosen here for convenience, to adequately suppress LFV. The choice 
for the elements of X is fixed, at the level of orders of magnitude, by the relation (39). The 
transformations needed to diagonalize the lepton mass matrix are given in eq. (38). 

Now we can complete our discussion concerning the Wilson coefficients that are bilinear 
in the spurions X, X. Coefficients {Nx, Ny) = (2, 2) of the vector operators are of the type 
XAX'^ or XAX"^ where A depends on c, c and Y^. From tables 2 and 3, we see that Y^ 
enters A only in the combinations Y^Y^^ or Y^^Y^. Therefore, up to an overall coefficient of 
order one, in the basis where X and X are diagonal and real we have 


When we move to the charged lepton mass basis, we get: 

The most stringent limits of table 7 arise from p“Au —)■ e“Au and require 


y 


^X^x, 
X 


2X^X, 


A" 

X 


~ y ~—~ < 5.0 X 10 TeV 


-2 


(46) 

(46) 


(47) 


which translate, respectively, into A > 2.2 TeV. Here and in the rest of this section A 

stands for either M or niQ. The limit on the decay t ^ gives rise to a similar bound: 


yi XrXi^ ~ y^ < 1-3 X 10 ^ TeV ^ , (48) 

resulting in A > ^ 1.9 TeV. Coefficients of the type {Nx,Ny) = (2,0) can arise for the 
contact operators of the type llqq and we expect, in the mass basis: 

i^lq’ {^tu,d)ij ~ XiXj , {Ceg)ij, {Ceu,d)ij ~ X^Xj . (49) 

The limits of table 8 from p“Au —>■ e“Au require 


X,X, 

A 2 


X^X, 

A 2 


yA^ 


< 1.0 X 10"® TeV 


7-2 


(60) 


17 








and we get A > h.Q/^Jy TeV, respectively. 

We are left with the coefficients that are quadrilinear in X, X. In the IFV scenario the 
off-diagonal elements of the dipole operator are dominated by the coefficients 

iXYiXlXX% , (C„)y = ^ (XXtXY^X%, , (51) 

which, in the lepton mass basis are approximately given by 

{Ce'f)ij ~ - XiXj , ^ XiXj . (52) 

The limit 

{Ce^)e^e ^ 2.5 x 10-^° TeV^ 

h? 

corresponds to a bound A > (0.6/y^) TeV. Similarly, from 

< 2.7 X 10-"^ TeV-' (54) 

A2 

we get A > (1.1/^/y) TeV. We have checked that all the other limits listed in table 5-8 do 
not lead to more restrictive bounds on A. For instance, the bound on the contact operators 

^errfj, 

< 1.2 X 10“^ TeV-' , (55) 

A2 

produces the bound A > {0.3/y) TeV. The bounds from the vector operators scale as 
while those from the dipole or contact operators scale as 1/y or 1/^/y. Therefore a new 
physics scale A around few TeV is still acceptable provided the Yukawa coupling y is close 
to one. 

A new physics scale A 1 TeV is too large to allow, in the present framework, for an 
explanation of the central value of the Aa^ anomaly. Indeed, in the lepton mass basis we 
have 

~ = 3.8 X 10-® (56) 

and to match the required value, eq. (32), we need A = 0.56 TeV. Since the contribution to 
AOju scales with the inverse square of A, choosing A = 1 TeV gives Ao^u = 9 x 10-^°, less 
than one third of the central value of the current anomaly. Concerning the electron EDM, we 
assume as in MFV that the sources of CP violation and LFV are the same. Since the Wilson 
coefficients of the dipole operators with Nx = 2 are aligned in flavour space with the mass 
operator, we identify in (51) the dominant coefficients that can contain non trivial phases. 
We estimate 

/™(C“) =. ^ , (67) 

which, for A > (3.1Q/^Jy) TeV, respects the bound of eq. (30). 

In summary, even in the limit of “wrong” Yukawa coupling negligibly small, an anarchic 
requires a scale of new physics well above 10 TeV. One way to lower this bound consists 


(53) 
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in mimicking the case of MFV, where there is a single non-universal spurion, either X or 
X. In this case, when neutrino masses are neglected, LFV is absent. By discussing a large 
class of possible LFV effects, we have shown that MFV is not the only possibility to reconcile 
LFV bounds with a new physics scale around few TeV. In the example analyzed here, both 
X and X are non-universal and represent potential sources of LFV, which is allowed also in 
the limit of vanishing neutrino masses. The corresponding bounds on A approach the TeV 
scale, provided the Yukawa couplings are close to one. 

4 Two-site model 

In the previous sections, we focused on PC scenarios for charged leptons from a general 
perspective, exploiting a spurionic analysis. Now, instead, we consider a specific simplified 
composite Higgs model, the so-called two-site model [37]. Its relevant features are: 

(i) The gauge group is = G®' x where 

= [SU{2)l X U{1)yY^ , = [SU{2)l x SU{2)r x • (58) 

The group for has been chosen in order to provide a custodial symmetry. In 

table 9 is summarized our notation for generators, boson fields and coupling constants 
associated with each simple subgroup. 

PC is assumed to arise from an unknown dynamical spontaneous symmetry breaking 
mechanism, that take place at higher energies and that breaks G®* x G®°™p into the 
diagonal group (which can be recognized as the electroweak gauge group). This spon¬ 
taneous symmetry breaking mechanism, that in turns triggers the PC scenario, will be 
effectively described by linear couplings between elementary and composite bosons. 



Subgroup 

Generator (s) 

Field(s) 

Coupling 

el 

SU{2)l 

T^L,ei^ a =1,2,3 

W“. el 

V 

gf 

U(l)y 

Y 


gf 

comp 

SU{2)l 

r“Ucomp^ a = 1,2, 3 

ppa, comp 

^comp 

SU{2)n X U{l)x 

c°mP+ yiv) 

Dcomp 

^comp 

comp ^2i^, comp 

T^1,2, comp 

comp 

9 2 

comp _ 

Dcomp 

comp 

g2 


Table 9: Gauge subgroups and their associated generators, boson fields and couplings. The 
normalization of the and generators has been chosen to match the S'O(IO) GUT 

normalization of the hypercharge, Ygut = y |^- 


(ii) The fermionic sector of the model includes three families of chiral fermions charged 
under G®* and three families of vector-like fermions charged under G®°™p. PC for 
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Elementary 

Composite 

SU{2)l 

U(l)y 

SU{2)l 

SU{2)r 

U{l)x 

('Ll 

2 

1 

2 

1 

1 

0 

^Ri 

1 

-1 

1 

1 

0 

Li 

1 

0 

2 

1 

1 [i 

2 ■ V 2 

Ei 

1 

0 

1 

1 

-1 • 

y 2 


1 

0 

2 

2 

0 


Table 10: Particle content and quantum numbers of the two-site minimal model. The index 
i = 1, 2,3 runs over three families for each representation. Lower case letters denote elemen¬ 
tary fields, capital letters denote composite fields. The ‘tilde’ apex denotes SU{2)l singlets, 
in order to distinguish them from the doublets. 


fermions is realized through linear mass-mixing terms between elementary and com¬ 
posite fermions. 

(iii) The Higgs sector consists of a real bidoublet {(f, (p) charged under [SU (2 )l x SU (2)ij]“™P, 
that will be identified with the composite Higgs field, interacting only with the com¬ 
posite fermions. 

Table 10 summarises the quantum numbers for fermions and Higgs doublet. 


We are ready now to introduce the Lagrangian of the two-site model: 

-^el T -^comp T -^mix ■) 


with 


1 ^ - 


1=1 


1 


= - t(p“ )"+ \DM‘ - n<p) 


a \2 
^ \r iiu J 

3 

+ ^ - mi)Li + - fhi)Ei 

1=1 
3 

- ^ (v^^jlRi^ELj + Y^-jLiiipERj^ + h.c. , 

hj = ^ 

C ■ - M (^A ^ ' 

“ 2 ^comp"^ + 2 \gComp^^ 

3 

^ ^ -|- h.c. , 


(59) 


(60) 


(61) 


(62) 
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where lower (upper) case letters denote elementary (composite) fields and the ‘tilde’ denotes 
SU(2)l singlets. Moreover, among the heavy vector bosons, we distinguish between those 
mixing with the SM gauge bosons, p* = {W*, S))}, and those that do not, 

In £ei we can recognize a SM-like Lagrangian (with collectively denoting the field 
strength tensors for the elementary gauge bosons) without the Higgs doublet, that is now 
a composite particle. Ccomp contains mass terms for composite leptons, kinetic terms for 
composite leptons and bosons (we have collectively denoted the field strength tensors with 
p“y), their interactions, the Higgs sector and Yukawa interactions. Finally, £mix contains 
linear mass-mixing terms among elementary and composite particles. We have also included 
in ^mix a mass terms for composite vector bosons. /Imix explicitly breaks G®' x down 

to the diagonal subgroup, and effectively reproduces the mechanism of PC. 

From the above Lagrangian we can easily recognize the flavour symmetry group of eq. (1) 
and convince ourselves that promoting Y^, Y^, m, fh, A and A to spurions with the trans¬ 
formation properties of eq. (3) we actually restore Gf. 

The Lagrangians of eqs. (60), (61), (62) are expressed in the elementary/composite basis. 
For the bosons, in order to switch to the mass basis (before EWSB), we diagonalize the 
mass-mixings in by the following field transformations: 


\ J V COS0 ) \ Pn ) 

Pn ^ Pn , 


tand = 



^comp ’ 


(63) 


At this stage the fields are massless, while p* and p^ have masses of order M*. Taking 
into account EWSB, the diagonalization of the mass terms becomes rather involved. To a 
good approximation, eq. (63) still diagonalizes the boson fields, while the mass matrices for 
the fermions read: 
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(64) 


With this notation, one can then perform a suitable rotation that brings to the mass basis: 







(65) 


Nh ^ Un N. 




where Vr^r, Ul,r are defined in order to have VIMsVr and diagonal. 

Using rotations (63), (65) one can recast the Lagrangian (59) in terms of lepton mass 
eigenstates. Using an approximate expression for the rotation matrices Vl,r, Ul,r, assuming 
universal masses for the heavy leptons, (m* = m,rhi = m), and retaining the leading terms 
relevant for FC neutral-current (FCNC), one gets: 

£fcnc = T (^(XY^)„euER, + (YG^EuiRi) + h.c. 
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XY/Y^^X^ 


yr- 


2m? 




mm 


-XY*Y^ 
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^Li^Rj + h.C. 




+ ^2/?“ +h.c. 

+ 1?-^“ (/(XyRYx'heuYeL, -/(XY^*YiX')Yf^YSR,) 

t : \ II i lli J 

- |(B; - B„ + w;-^) (XX%euYeL, + S.(B; - B,) (XX%,f.R,YeR, 

+ |(B; - B, + W;-’>) X.YuYEl, + g.(B; - 4) YlaYEn, + li.c. (66) 


From Lagrangian (66) one can easily read some of the most relevant features of FCNC in 
this class of models: 


• The first line of eq. ( 66 ) accounts for FCNC interactions among fight- and heavy-leptons 
and the Higgs. They contribute to low-energy processes such as £j —?• ij'j at loop-level. 

• The second line of eq. ( 66 ) refers to FCNC interactions among light-leptons and the 

Higgs. As we know from previous sections, FCNC effects can arise in this case through 
the dim -6 operator {Qe/ij = /)/Li&Rj/) and this explains the appearance of the 

factors jm? ^ I'm? and jm'fh after EW symmetry breaking. Such Yukawa inter¬ 
actions generate FCNC Higgs decays like (p —>■ r/r, tree level contributions to processes 
sensitive to four fermion operators such as p —>■ 3e and \xN —>■ eX and loop induced 
effects to A ^ 

• The third line of eq. ( 66 ) describes the interactions among light- and heavy-leptons and 

the Z boson. Since the operators eu$ELj and Em/eRj are not SU{2)l invariant, they 
can be generated only after the EW symmetry breaking through the SU ( 2 )^ invariant 
dim-5 operators lutERjLp and Lm/eRjif. This explains the factors of v/'fh and v/m. 
These interactions will contribute to at loop-level. 

• The fourth line of eq. ( 66 ) contains ECNC interactions among light-leptons and the 

Z bosons. As discussed in previous sections, such effects can arise by means of dim -6 
operators bilinear in the Higgs doublet like (Qtpe)ij = /){?Ril^eRj). As a result, 

after EW symmetry breaking, we generate the operators eu/^Lj and eRi/cRj which 
are suppressed by /m? and I'm? factors, respectively. The leading effects induced 
by these terms are the tree level FCNC decay modes Z E/j as well as —?• 3£j and 

jj,N —>■ eN. 

• The fifth line of eq. ( 66 ) describes SU{2)r invariant interactions between heavy gauge 

bosons and light-fermions. Also in this case we can induce FCNC processes at tree level 
such as C Mj and jiN —?• eN and loop induced effects to —?• ^^ 7 . 

• The sixth line of eq. ( 66 ) refers to FC interactions among heavy gauge bosons, heavy- 
and light-leptons. Effects to low-energy observables are induced by the loop exchange 
of heavy gauge bosons and leptons. If the heavy leptons of different generations were 
degenerate, ECNC effects would vanish according to the GIM-mechanism. However, 
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the GIM cancellation is broken by non-universal mass splittings of order X'^X and 
XX. The latter point has been overlooked in the literature so far. 

Concerning the flavor structure of the various interaction terms, we remember that each Higgs 
electroweak doublet ip or its vacuum expectation value v occur accompanied by a power of 
or Y^, while the lepton fields and cr come with X and X, respectively. 

Before discussing LFV in this model we comment on the issues of renormalizability, gauge- 
invariance and UV sensitivity. The two-site model is a non-renormalizable effective theory 
and therefore it is valid only up to energies of the order of an UV cut-off Auv Such a cut-off 
has been estimated in ref. [37], where it has been found that Auv ~ SnM*/g*, by analysing 
self-interactions in the composite sector. The non-renormalizability is a consequence of the 
explicit breaking of the gauge symmetry x by the heavy vector boson masses, 

by the Yukawa couplings Y^ and Y^ in the composite sector and by the mixing between 
elementary and composite fermions A and A. 

This explicit breaking raises the question of the reliability and consistency of our results. 
We can always promote our effective Lagrangian to a gauge-invariant theory by interpreting 
the sources of gauge symmetry breaking as spurions. Providing the spurions with suitable 
transformation properties under the gauge group we can recover the gauge invariance under 
the full local group G®' x G‘^°™p. Once we treat the spurions as dynamical fields, this procedure 
defines a (non-unique) embedding of the effective theory in a possible UV completion. The 
spurions include two types of degrees of freedom: the would-be Goldstone bosons, eaten by 
the heavy gauge vector bosons of the composite sector through the Higgs mechanism, and 
physical scalars By working in the unitary gauge, which we adopt in our computations, a 
generic amplitude comprises two separate contributions: one coming only from the exchange 
of the physical polarizations of the gauge vector bosons and one including the exchange 
of some physical scalar degrees of freedom $*. Clearly our computation retains only the 
first one, while the second is missing and makes our results sensitive to the details of the 
UV completion. If the masses of the extra scalars are close the cut-off scale Auv, we 
expect that the contributions we are neglecting in our estimates of the Wilson coefficients 
are generically of order [M/AuvY and/or (m,7h/ AuvY ■ 

The embedding of the model in a UV completion also shows that it is not consistent to deal 
with the HF case, within the effective theory. Indeed to keep the masses of the elementary 
leptons non-vanishing in the HF limit, we have to consider at the same time large (A, A), 
such that the ratios {X,X) remain constant. In a UV completion A, A are proportional to 
VEVs that break the gauge symmetry G®* x G'^°™p and contribute to the masses of the heavy 
gauge vector bosons of the composite sector. Therefore, barring tuning of the parameters, 
we cannot make m,m S> M. 


5 Lepton flavour violation in the two-site model 

In this section we will present our results for the Wilson coefficients of the various LFV 
operators considered in sec. 2, in the context of the two-site model introduced above. We 
work at the leading order in the loop expansion and we pay particular attention to the 
spurionic structure of the coefficients. Throughout this section we assume universal masses 
for the heavy leptons: m, = m and m* = fh. 
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Figure 1: Feynman diagrams for LO contributions to £j ij'y with loop-exchange of SM 
bosons h, Z, W and heavy fermions (left) and heavy gauge bosons and heavy fermions (right). 


5.1 Dipole operators 

Dipole operators are among the most important operators in the leptonic sector as they 
generate LFV radiative decays such as fi ^ ey, leptonic EDMs and (^r —2)s. We will consider 
the one-loop contribution to dipole operators in the two-site model, assuming the validity of 
the perturbative expansion. The result of our computation of the full one-loop contribution 
to the electromagnetic dipole operator is shown in the appendix B. Given the fact that the 
composite Higgs scenario is naturally characterised by a strongly interacting regime, the use 
of perturbation theory can be questioned. Here we will restrict ourselves to the portion of 
the parameter space that is compatible with perturbation theory. This allows us to consider 
moderately large coupling constants of the composite sector, with loop factors still providing 
a sizeable suppression. The same considerations apply to processes receiving non-vanishing 
contributions already at the tree level. Within the two-site model, dipole amplitudes receive 
leading contributions from the virtual exchange of i) light SM bosons (h, Z, W) and heavy 
fermions and ii) heavy gauge bosons and heavy fermions. Examples of diagrams for the 
above classes of contributions are shown in fig. 1. In particular, summing over the h, Z and 
W amplitudes, we find that the dominant effects for class i) are given by 


h? Gdyr^ mm R l r 

in agreement with the result of refs. [35,62]. As we can see, the amplitude of eq. (67) has 
precisely the same spurionic structure as the one of eq. (6). In particular, according to our 
spurionic classification of sec. 2, it turns out that order (Xy,Xx) = (3,2), 

where Ny and Nx are the orders of the expansion in (Y^,Yl) and (X, X), respectively. 
Interestingly, we observe that contributions of order {Ny,Nx) = (3,2) containing only the 
Yukawa Y^, which would be allowed by the flavor symmetries of our model, are absent, to 
one loop order. 

The leading effects for class ii) start from the order (Ny, Nx) = (1,2) and read 




eiJB+B* 
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2567r2M2 L 

e gl 

2567r2 M2 


XY’^Pf^(v, z) + XYlX'g(y, z) 

xY,;x^fr(i/), 
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( 68 ) 

(69) 


















where M stands for a common heavy boson mass, y = z = and the loop 

functions are defined in the appendix. Notice that the first term of (Cg-y)and are 

aligned with the mass operator as long as m,rh ocl, at least. Therefore, as already pointed 
out in the literature [12], they cannot induce neither flavor nor CP violating effects after 
switching to the mass basis for the SM leptons. Yet, we stress that they contribute to the 
leptonic g — 2. However, we point out here that the second term of {Cg ^)not discussed in 
the literature to our knowledge, is not aligned with the mass operator and therefore generates 
flavor and CP violating effects. 

From eq. (67) one can easily find that, in the anarchic scenario, the bound from BR(/^ —)■ 
ey) imposes that -s/mm/(Y) > 10 TeV. In order to relax such a strong bound while keeping 
anarchic, we analyze here in great detail the solution with = 0, as already done in the 
model-independent analysis of section 2. 

Setting Y£ = 0, the NLO effects stemming from class i) are given by: 


(cWt 

A2 


2567r2 


XY^X^XY^^Y^X^ 






-R 


A2 


C^YY^Y^^X^YY^Y^ + ^^^YY^Y^YY^^Y^Y^ 


2567r2 
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{m? — fh?) 


rYY^Y^V^Y^^Y^Y^ + 8 YY^Y^V^Y^^^ 


(C„) 


w 




256^= 


(70) 


(71) 

(72) 


where x = rX/fX. The first (second) terms of {Cg^)^,{Cg^)^ and (C'g.y)^ arise from dia¬ 
grams such as those shown in the left (right) plot of fig. 2. These contributions are of order 
(Ny,Nx) = (3,4) and occur through a chirality flip implemented in the internal heavy- 
fermion line (left plot) and external light-fermion line (right plot). On the other hand, the 
last term of {Cg^)^ as well as the third terms of {Cg^)^ stem from dim-8 operators such as 

^(?iV3)(a . F)X Yi X<eK , (73) 

and the corresponding diagrams are shown in the left plot of fig. 3. Such contributions are 
of order (Yy, Nx) = (5, 2) and turn out to be suppressed by a factor of /fh? compared to 
dim-6 contributions. As a result, the parametric ratio between dim-6 and dim-8 operators is 
of order X'^fh?/Y^v^ showing that both kind of operators might provide the dominant effects 
depending on the model parameters. Let us stress that contributions of order (Yy, Yx) = 
(5,2) can arise also at dim-6 level if the extra Higgses of eq. (73) close in a loop instead 
of getting a vacuum expectation value. Therefore, we expect that one-loop induced dim- 
8 operators with (Yy,Yx) = (5,2) dominate over two-loop induced dim-6 operators with 
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Figure 2: NLO Feynman diagrams for £j —?• ij'y of order (Ny, Nx) = (3,4) with loop-exchange 
of h,Z,W and heavy fermions. 



Figure 3: NLO Feynman diagrams for £j^ arising from dim — 8 operators. The diagram 
on the left is of order (Ny, Nx) = (5, 2), the one on the right of order (Ny, Nx) = (3, 2). 


{Ny,Nx) = (5,2) provided j'ffZ' > I/IGtt^. The last term of (C'e.y)^ stems from dim-8 
operators of the form 


■ F)X X'Ir . 


IGtt^ 


(74) 


and the relevant diagram is shown in the right plot of fig. 3. Finally, we observe that there 
are not effects for {Ce-y)^, from dim-8 operators since our model doesn’t include heavy right- 
handed neutrinos. 

We discuss now NLO effects arising from class ii) which are given by: 
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(76) 


The first two terms of and {Ce'y)y^* are of order {Ny,Nx) = (1,4) and arise from 

the one loop-exchange of heavy fermions and bosons. The relevant Feynman diagrams are 
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Figure 4: NLO Feynman diagrams for —)■ f'j-y of order (AV, Nx) = (1,4) with loop-exchange 
of heavy gauge bosons and heavy fermions. 


shown in fig. 4. On the other hand, the last contributions to and are of 

order {Ny, Nx) = (3, 2) and stem from dim-8 operators such as, for example, 

■mr^(mtm)“‘v'yyS,VeK . (77) 

The relevant Feynman diagrams are shown in the right plot of fig. 3. 


5.2 Scalar operator 

The Yukawa interactions for charged leptons are modified after we integrate out at tree 
level the heavy fermions. Setting = 0, we find that: 


/j _ 

-Ly = -^eLyeeR + h.c. , 

where in the mass basis for the charged leptons we have 

» = »f' - ^ XY^YilX'yr - ^ yfXYi^Y^P , 


(78) 


(79) 


with 

yr = XY*X^ . (80) 

Therefore, the coefficient of the dim-6 operator = C'e<^((/?V)(/Leij(p) reads 
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Notice that for = 0 the corrections to yf^ are proportional to yf‘ 
where the Higgs is a pseudo-Goldstone boson [63-65] . 


(81) 

itself, as in scenarios 


5.3 Vector operators 

The Z boson interactions with charged leptons are also modified after we integrate out 
at tree level the heavy fermions. In particular, we find that 

= — e (ghPi + gnPR) 1 (82) 
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where ql and qr are defined as follow 


= -l + si, + ^XYiYi'XK 




Switching to the coefficients C*® and C^e we find that: 


- -XYXYX^X^ , 

2m2 -R ^ 


5.4 Contact operators 


cS’ = 0, 


A2 2m2 ^ 


Contact operators induced by the exchange of heavy gauge bosons can be easily derived 
starting from the coefficients and entering the interaction Lagrangian of heavy 
gauge bosons and SM leptons which we report in the appendix. Integrating out the heavy 
gauge bosons at tree level, we find that 
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where higher-order effects suppressed by additional factors of [gsM/g*Y, v'^/m?{7v?) and 
X2(X2) have been neglected. 


6 Phenomenological analysis 

In this section, we evaluate the most relevant low-energy processes in the charged lepton 
sector in the context of the two-site model by making use of the Wilson coefficients derived in 
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the previous section. One of the most interesting features of the two-site model is the absence 
of the Wilson coefficients {Ny, Nx) = (3, 2 ) from dimension six dipole operators, at least in 
the one-loop approximation. Such coefficients are particularly dangerous for LFV. As we saw 
in section 3, their presence is not compatible with anarchic Yukawas if the compositeness 
scale is close to 1 TeV. We do not know if this feature of the model is an accident of the 
one-loop approximation or if it persists at higher loop orders. In what follows we will neglect 
the (Ny,Nx) = ( 3 , 2 ) contribution to dimension six dipole operators, assuming that higher 
loops provide a sufficient suppression. Our purpose is to check if under these conditions 
anarchic Yukawas are still viable or not when we consider a compositeness scale around 1 
TeV. 

Among the most interesting LFV channels are fi —?• ey, /r —?• 3 e, /i ^ e conversion 
in Nuclei as well as r LFV processes. However, hereafter, we focus on processes with an 
underlying /r —?• e transition since they are the best probes of composite Higgs models with 
anarchic Y^. Concerning flavour conserving processes, we are interested in the electron EDM 
and the muon g — 2. The current status and future experimental sensitivities for the above 
processes are collected in table 4. As recalled in section 3, there is a ~ 3.5(J discrepancy 
between the SM prediction and the experimental value of the muon g — 2. 


6.1 £-1 —y 

The dipole transition ^^7 is responsible for both LFV radiative decays (when i ^ j) 
like /i — 7 - ey and flavor conserving processes like the electron EDM and the muon g — 2 when 
i = j = e ox ji, respectively. The branching ratio for the process /r ^ ey can be written as 


BR(/i ey) 


247r2 



(87) 


Radiative LEV transitions are tightly related to the magnetic and electric leptonic dipole 
moments which are also extremely sensitive probes of new physics. In particular, one can 
find that 

4 = - ^ Im (C“) , Re (C)-;) . (88) 

In concrete scenarios as our two-site model, Ao^, and BR(£ —>■ Ay) are expected to be 
correlated. However, such correlations depend on the flavor and CP structure of the couplings 
which are unknown. In our discussion, we assume order one CP-violating phases and anarchic 

Here we provide the dominant contribution to /r —>■ ey which arises from dim -8 operators. 
Hocusing on the HB scenario (M ^ m, m) with anarchic Y^ and Y = Y, it turns out that 


BR(;U —)■ ey) 


647r {w? — ’ 


(89) 


where we made use of eq. (40) to eliminate Y, Y and where now stands for an average 
element of the anarchic matrix Y^^y Notice that the above expression is valid only in the 
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Figure 5: Branching ratio of ^ 07 as a function of the heavy fermion mass m. 


mass insertion approximation which requires that |m — m| ^ As an example, if 

|m^ — m^l one can find 


BR(// 07 ) 3 X 10 


/ 1.5TeV \ 
\ m J 


I V* 1^ 
\^r\ ^ 


(90) 


implying that p ^ ey saturates its current experimental bound for m ~ 1.5 TeV and 
~ 1. This is conhrmed by our numerical results shown in hg. 5 where we have reported 
the predictions for BR(/i ^ ey) as a function of the heavy fermion mass m. In fig. 5, as well 
as in all other plots, we have assumed anarchic and Y£ and X = X so that the relevant 
flavor mixing angles entering /r —?• e transitions are Xg/X^ = ^ The lower 

(upper) red line in the left plot refers to the case where Y^ = Y^' = 1 (Y^ = Y^ = 2 ) while the 
lower (upper) black line corresponds to Y^ = 0 and Y^ = 1 (Y^ = 2 ). The most prominent 
features emerging by this plot are: i) in the case of Y^ = 0 {¥[ ^ 0 ), the bound on /i —?• ey 
is satisfied for m > 1 TeV (m > 10 TeV), ii) for Y^ = 0, BR(/i ^ ey) ~ while for 

7 ^ 0, BR(/i ey) ~ and this explains the different growth with Y^,Yl and 

the decoupling properties with m of the various curves. In the right plot of fig. 5 we show 
the branching ratio of p —>■ ey as a function of the heavy fermion mass m for Y^ = 0 and 
various values of Y^. 

A quite similar behavior is expected for the electron EDM. Indeed, if Y^ 7 ^ 0 and assuming 
0(1) CP-vioIating phases, it turns out that 
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(91) 


and therefore the electron EDM poses a severe constraint on the heavy fermion scale at the 
same level of fi ^ ey. Setting Y£ = 0, the EDM bound is satisfied with much lighter masses. 
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Figure 6: Electron EDM dg as a function of the heavy fermion mass m. 


In particular, we find that 
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(92) 


where in the last equality we have assumed that \m? — m^| ~ m?. 

In fig. 6 we show our numerical results for dg. Indeed, fig. 6 represents the analogous of 
fig. 5 for the electron EDM case and similar conclusions drawn for p ^ eq apply here too. 


6.2 fjL —^ 36 

The process /r ^ 3e receives more contributions than fi ^ ej as it is sensitive also to 
four-fermion operators. In particular, starting from the general expression of ref. [53,66] and 
neglecting scalar and non-dipole photonic contributions which are negligible in our model, 
we find the following expression for its branching ratio 


BR(/i —7- 3e) = 
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8ve ^ 
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^ {‘^CvLL + Cvlr) C'e’T* + {‘^CyRR + CvRl) 


m 


+ ^ - ^ ) (ICICI" + ICICI") 


mt 


where 


( 93 ) 


CvLL = (24- - ^cS"" + qr. CvRR = 24CJC + Cjr. 
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CvLR = + cr' . CvRL = (2sl. - 1)0* + Cjr" , (94) 

and the explicit expression for the he Wilson coefficients Cab can be found in sec. 5. In 
the above equations, we have neglected scalar four-fermion operators since they are very 
suppressed in our model. We remind that, whenever the dipole operator is dominant in 
fj, ^ 3e, there exists a model independent correlation between the branching ratio of /r ^ 3 e 
and jjL ^ given by 


BR(/i — 7 - 3e) 
BR(/i — 7 - ey) 


(^el 

Sir 
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m,, 
log—I 



6.6 X 10 ^ . 


(95) 


As a result, the current MEG bound BR(p —>■ ey) < 5.7 x 10“^^ already implies that 
BR(p —>■ 3e) < 3.8 x 10“^^. However, as shown in our numerical analysis and illustrated 
in fig. 8 , p ^ 3e turns out to be dominated by non-dipole operators and the correlation 
of eq. (95) is significantly violated. This is a relevant result, as within composite Higgs 
models with 7 ^ 0, as well as in supersymmetric scenarios, eq. (95) holds to an excellent 
approximation. Therefore, the two-site model with = Q can be disentangled among other 
models if both /r ^ 3e and /r ^ ey will be observed. 

In particular, in the HB scenario (M ^ m, in) with anarchic Y^, m = rh and X = X, we 
find 


BR(/i ^ 3e) 
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(96) 


implying that fx ^ 3e saturates its current experimental bound for m ~ 1 TeV and Y^ ~ 1. 
This result is well reproduced by our numerical analysis shown in the left plot of fig. 7 where 
the behavior of BR(/i —)• 3e) as a function of the heavy fermion mass m is shown for different 
values of Y^. Notice that BR(/r ^ 3e) ~ /rri^ and therefore it has a much milder growth 

Y^/mP' (in the case 


with Y^ and a slower decoupling with m compared to BR(p ey) 


of 12 = 0). 


6.3 ijlN —^ gN 

As the last process governed by an underlying [x ^ e transition we consider /r — e con¬ 
version in nuclei. As in the /x ^ 3e case, this process receives effects from both dipole and 
four-fermion operators. In particular, the fx — e conversion branching ratio is defined as 

BR{ix-N e-N) = , (97) 

^capt 

where cucapt is the muon capture rate while tUconv is given by the expression 



32 











-10 




Figure 7: Branching ratio of 3e (left) and fi Au e Au (right) as a function of the 
heavy fermion mass m. 


which we have obtained starting from the general expression of ref. [56] and setting to zero 
the scalar and non-dipole operator contributions which are negligible in our model. The 
coupling constants Al, Ar, and g^y^ in eq. (98) are defined as 
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The flavour-conserving interactions of the heavy vector resonances with the light quarks 
have been derived from the respective couplings with the light leptons simply rescaling them 
according to the different charges of the quarks under Y and T 3 . Finally, the quantities D and 
VP,n refer to the overlap integrals between the wave functions and the nucleon densities [56] 
which we report for few relevant nuclei in table 11. An inspection of eqs (97),(98) and table 11 
shows that at present g~Au e~ Au is the most sensitive probe of new physics among the 
various ji — e conversion in nuclei processes. In scenarios with dipole dominance, the following 
model-independent relation holds 


BR(;U Au —>■ e Au) D 
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7 ^_ 7 __ e' 

BR(yU ^ eq) 1927r2 ujcapt 


■ 3.8 X 10 


-3 


(106) 
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Nucleus 

D 

y(p) 


^capt(lOS-l) 

Au 

0.189 

0.0974 

0.146 

13.07 

Ti 

0.0870 

0.0399 

0.0495 

2.59 

A1 

0.0362 

0.0161 

0.0173 

0.7054 


Table 11: Overlap integrals from ref. [56]. 


and therefore Au —>■ e~Au) < 2.2 x 10“^^ after imposing the bound BR(/r —>■ 67) < 

5.7X10“^^. However, as in the case of /r —>■ 3e, in turns out that in our model BR{fi~Au —>■ e~Au) 
is dominated by non-dipole operators and the correlation of eq. (105) is not at work. This 
is shown in fig. 8. In particular, in the HB scenario (M ^ m,m) with anarchic Y^, m = m 
and X = X, we find 


BR(/r Au e Au) ~ 


^capt 


^ rraerra^ ^ 


|p^O)|2 


4 X 10“^^ 



(106) 


and therefore the current experimental bound is saturated for m ~ 3 TeV and ~ 1. 

These expectations are fully confirmed numerically as shown by the right plot of fig. 7 
where we we report BIl{fi~Au —?• e~Au) as a function of the heavy fermion mass m for 
different values of Y^. Since both BB{ii~ Au —>■ e~ Au) and BR(/r —>■ 3e) are proportional to 
/rrA, the consideration done for the case of /U —>■ 3e hold here too. However, we point 
out that in our model BB{ii~Au —>■ e~Au) is much enhanced with respect to BR(/r 3e) as 
opposite to scenarios with dipole-dominance where BR(fi~Au e~Au) ~ 0.6 x BR(/i —?• 3e). 

As a result, the simultaneous observation of these processes would enable us to disentangle 
among the underlying theory at work. 


7 Conclusion 

We have reanalysed the bounds on the compositeness scale derived from LFV and CP 
violation in the lepton sector, in the framework of PC. In the generic case of anarchic Yukawa 
couplings, such bounds are known to be quite severe, stronger than the analogous bounds 
from the quark sector and requiring a compositeness scale above 30 TeV. In this work we 
focussed on the case of vanishing neutrino masses and vanishing “wrong” Yukawa couplings, 
in the hope of minimising FCNC and CP violating effects. We performed a general effective 
operator analysis, where the Wilson coefficients of the relevant dimension-six operators are 
determined by a flavour symmetry group and by a set of spurious. We have considered 
all lowest dimensional operators leading to LFV violation and CP violation that can be 
constructed with the Higgs doublet and the SM leptons. We have formally expanded the 
Wilson coefficients in powers of the Yukawa couplings of the composite sector and in powers 
of the elementary-composite mixing terms. Our expansion includes, for the first time, terms 
quadrilinear in the elementary-composite mixing (Y, X ). By exploiting the known limits on 
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Figure 8: Branching ratio of /r —>■ 3e (left) and [i Au e Au (right) versus the branching 
ratio of 67 for = 0. The case of dominance of the dipole operator is shown in yellow. 

the Wilson coefficients we have shown that, even in the case of vanishing “wrong” Yukawa 
couplings, the anarchic scenario is not compatible with a compositeness scale of 1 TeV, barring 
accidental cancellations not incorporated in our general spurion analysis. We have also 
shown that there are schemes where a TeV scale can be accommodated, without necessarily 
reproducing the case of MFV. In the example we have provided this is achieved in a semi- 
perturbative regime, where the universal Yukawa coupling y of the composite sector is close 
to unity. In this respect it is interesting to note that while the bound on A^vp coming from 
operators bilinear in (X,X) scales with y/y, the limit from the operators quadrilinear in 
{X,X) scales with 1 /y, thus making essentially impossible to lower Ajvp below the TeV. 
This also shows the importance of keeping terms quadrilinear in the mixing in the expansion 
of the Wilson coefficients. 

We have also derived the low-energy effective Lagrangian relevant to LFV and CP vi¬ 
olation in the lepton sector in the two-site model, that includes heavy vector-like fermions 
as well as heavy spin-one particles. We focused on the case = 0 , finding a strong sup¬ 
pression of LFV effects compared to a generic composite Higgs scenario where Y^ ^ 0 . In a 
perturbative analysis where we evaluated the relevant amplitudes at the LO, we found that 
the /i —7- 67 transition is dominated by dimension 8 operators, at variance with our spurion 
analysis where the most important contribution comes from operators of dimension 6. As a 
result the best probes of our scenario are /r —>■ e conversion in nuclei and the electron EDM. 
A1 least in that portion of the parameter space where perturbation theory is applicable, LFV 
allows a compositeness scale close to the TeV, even in the case of anarchic Yukawas. In 
this regime there are interesting relations among the various LFV transitions, which allow to 
disentangle the model from other possibilities. 
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A Contact operators 

In the following, we provide the full results for the coefficients and entering the 
interaction Lagrangian of heavy gange bosons and SM leptons 

r = ft. <V (cgPi + c^Pr) i , (107) 


MB*) _ ^ 
^Ll 2 


with pp = {B*, B,W^}. For the coefficient, we hnd: 

tan(d2)^ - (tan(02)^ + l)XX^ + (tan(d2)^ + 1)XX^XX^ 

f + 2) VA'YJi'jlOft - i(tan(02)2 + + h-c.) , 

: 777 / ^ 777 / 777 / / 


1 

i' 




XX<{ 1 -XX<) + ^^XY^Y, 


XY*Y*^ X^+h.c. 


2 cos 02 \ '' ^ fh? ^ ^ 2 mm 

cm) ^ tan(ffi)2 + (tan(ffi)2 + 1)XX^ - (tan(ffi)2 + 1)XX^XX^ 


2 taii(«i)Wxi-ryAl + htan(«i)" + 1) 

4 777/^ Z 


777 / 777 / 


XY^Yl^'X' + h.c. 


(108) 

(109) 


( 110 ) 


For the coefficient, we hnd: 


^(s*) 

^RR 


r<^B) 

^RR 


^RR 


/(tan(02)^ - (tan(d2)^ + l)XX^ + (tan(02)^ + l)XX^XX^ 

-i(2tan(«2)2 + i)YxyMM+ l)-GAKpyR-At+h.c.). (Ill) 

8 2 777 / 777 / / 

(XX^-XX^XX^ + -^XYX^YXX^+ ( -^YWV;^Y^+h.c.l), (112) 

cos02\ 4m2 ^ ^ \2mm ^ ^ JJ 

- (113) 
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Mb j Cj 

(1) (2) (3) 

Figure 9: One-loop diagrams contributing to the dipole operator in our model. In diagrams 
(1) and (2) we have ^ = 1,2, 3, in diagram (3) we have 5 = 1,2. 




Figure 10: Feynman rules needed for the computations. Here A,C = 1,2,3 and B = 1,2, 
while p = {Z, B*, B, W*’^}. The expressions for T, can be found in the appendix. 

B One loop contributions to the dipole operator 

In what follows, all lepton fields are meant to be in mass basis. We introduce a compact 
notation to collectively denote such fields, using an upper case latin index to distinguish 
among them: 


(■A/l)b = {j^L, Nl} , 


{£r)a — Er, , 
AfR = Nr , 


A = 1,2,3, (114) 

5 = 1,2. (115) 


The flavor index, when necessary, will be explicited with a lower case latin letter. For 
most of the expressions however this index will be omitted to simplify the notation. In those 
cases one should interpret the expressions as 3x3 matrices in flavor space. 

In figure 9 we see the different kinds of loop diagrams contributing to From the 

Lagrangian (59), after performing the rotations (63), (64) to the mass basis, one can read the 
Feynman rules of interest for these computations. We summarise them in figure (10), where 
the expressions for T, are: 




1 

o 

o 

_1 


1- 

o 

o 

o 

1_ 
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1_ 


_ 0 1 0 _ 


4?, 




ul 


1 0 0 
0 1 0 


Vr 


(116) 


37 














■ 0 

0 


_ 0 
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0 cot 6*1 0 


V, 


R 


= ul 


— tan 6*1 0 0 

0 cot 6*1 0 


_ T/t 

Cm 

■ 0 
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■ 0 
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2cw 
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14, 
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2cw - 



tan^ 62 



1 tan^ 62 
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1 

2 

14, 
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2 
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, (117) 

1 "l, 

(118) 

(119) 


Vr 


A 


= gi sec 6*2 l4 


W*’^ _ ^ T/t 

R - 2 


■ 0 

1 

2 

14, 

Af = g* sec 02 Vf 

■ 0 

1 

2 

1 



1 


■ 0 


A/'” = %Vl 

tan^ 01 

-1 

14 , 

-1 

0 


Z 

0 _ 


Vl, 




r = vl 


0 0 


0 


0 0 y, 


R 


0 y;i 0 




R 


( 120 ) 

( 121 ) 

( 122 ) 


where (yf* = gcotOi and cot 6*2. Remember that this is a compact notation in which 

the flavor indices are omitted. Then each element of the above matrices is really a 3x3 matrix 
in flavor space (such as Y^,Yl). 

We can now summarise the leading results for the different loop diagrams. In these 

2 2 

expressions, terms of order ^ or ^ (where A is either an heavy lepton or boson mass) have 
been neglected. The final results are: 
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(123) 
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16^/27r^v Ml 


2A^'iim^^A^ii + A^\2"^/f (y)A^2i + 


i^31 


+ mf^ (^-gA^n^Lll + ^il2^f2iy)^L21 + ^in'mf2{^WL3i 
+ + A£\2"^/2^(?/)Al2i + A^i3m/2^(^)A£'3i^ ^ , (126) 


gl 


32v^7r2w Ml 




21 


5M / -v=^t 


g^Lll^Lll 


Sl'i2>™/f(!/)S 


L21 


mf' 


gSi*„Ei„ + Ey,2n>/f(!/)S 


i ?21 


(127) 


J IJ 


where y = ■^, ^ = ■^, while in eq. (126) p' = {5*, B, W*’^} and the expression is the same 

p' p' 

for the three heavy bosons. Notice that the result is expressed in a block matrix notation. 
For example, taking eq. (123) and recalling that T is a 9x9 matrix, T 12 should be read as 
the 3x3 block of T in position ( 1 , 2 ). In this notation, the expression (Ti 2 h'^ 2 i)ij stands for 
the matrix element: Yfk=ii^ 12)kj- 

The expressions for the loop functions used in eq. (126), (127) are: 


/f(^) 

f 2 {x) 

fsix) 

fFi^) 


4: — 3x — + 6x log X 

(1 - xf ’ 
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10 — 43x + 78a:^ — 49a:^ + 4a:^ + 18x^ log x 
12(1 -a:)4 


(128) 
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C Loop functions 

The loop functions for dipoles at the LO read: 

fFiy , - + 


g 2{y + z){-A +3y + - Gylogy) 4 ( 2 ; - 2y)(-4 + 82 : + 2 :^ - d^^log^) 


fi{y,z) = 


_^/y^ 
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(1 - zf{y - z) 


-{y^z) 


friy) = 


A-27y + 2Ay‘^ - y^ - 6?/(l + 2y) log y 


{1-yr 

The loop functions for dipoles at the next-to-leading order read: 
4-7a: 
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